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On Abstract 

(— I This review summarizes recent experimental and theoretical progress in determining neutrino 

^5 mixing angles and masses through neutrino oscillations. We describe the basic physics of oscillation 

phenomena in vacuum and matter, as well as the status of solar, reactor, atmospheric, and acceler- 
ator neutrino experiments that probe these phenomena. The results from current global analyses 
of neutrino parameters are given. Future efforts that may improve the precision with which these 
parameters are known or probe new aspects of the neutrino mixing matrix are described. 



> 

P 1 Introduction 

(N 

The Standard Model (SM) of particle physics provides a renormalizable theory where charged fermions 
acquire Dirac masses through the Higgs mechanism. In this model the neutrinos stand out, lacking 
^ charges or any other additively conserved quantum numbers to distinguish particle and antiparticle, 
T-H and consequently having the possibility of both Dirac and Majorana masses, yet acquiring neither in 
> the SM. The structure of the SM - the absence of right-handed fields necessary for Dirac masses and 
^ of dimensional effective couplings necessary for Majorana masses - eliminates neutrino mass terms. 
^ Consequently the discovery 15 years ago of massive neutrinos supports the modern view of the SM as 
an effective theory, where for example Majorana neutrino masses would arise from the SM's unique 
dimension-five operator, and thus can be considered the most natural correction to that model. Hence 
neutrino masses represent an important opportunity to look for extensions of the SM and to determine 
the scales of the new physics that we now need to add to that model. Experiment has confirmed 
that there is a hierarchy of neutrino masses, just as there are differences between electron, muon, and 
tauon masses, though the overall scale of neutrino masses is radically different from that of the charged 
fermions, below ~ an eV rather than > MeV. This has implications for the scale of the new physics. It 
has also confirmed that there is fiavor violation among the neutrinos - the mass and weak interaction 
eigenstates are different - with the mixing angles being much larger than those found among the quarks. 
This has provided us with a powerful tool for probing neutrino properties, oscillations that couple the 
fiavors with significant amplitudes. 

The discovery of neutrino masses and mixing, as well as the changes in the mixing that accompanies 
neutrino propagation in matter, has stimulated a great deal of new experimental and theoretical activity. 
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Neutrino oscillation experiments have provided a reasonably clear picture of neutrino mixing angles and 
neutrino mass differences, though their are important issues not yet unresolved. Some of the remaining 
questions can be addressed with neutrino oscillations - e.g., the mass hierarchy, the Dirac CP phase, 
and the possible existence of additional neutrino generations - and some questions - the absolute 
scale of neutrino mass, the nature of that mass (its Dirac/Majorana character), and the Majorana CP 
phases - must be addressed by other techniques. This review summarizes the current status of neutrino 
oscillation studies and prospects for further exploiting this probe to test neutrino properties. 



2 Neutrino mixing and Oscillations 

The neutrino mixing matrix defines the relationship between the mass eigenstate and weak (flavor) 
eigenstate bases, 

I '^flavor ) T I l^mass ) • ( 1 ) 

We adopt the following parametrization of the neutrino mixing matrix: 
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where Cij = cosOij, Sij = smOij, and 5cp is the CP-violating phase. In writing Eq. ^ we have omitted 
the two Majorana phases because these play no role in neutrino oscillations. Note that the association of 
the CP-violating phase with 5*13 is simply a parametrization choice. In fact the CP-violating difference 
of the neutrino survival probabilities 

P{u^^lye)-P{l^^.^l^e), (3) 

is proportional to a product of the sines of all three mixing angles. 

Consider a particular flavor eigenstate, |z/j) = I]jTjj|i/j), created at some point in time by weak 
interactions, where the indices / and i denote flavor and mass eigenstates respectively. The subsequent 
propagation of this state is controlled by the free Lagrangian and consequently by the neutrino masses. 
As the neutrino travels downstream from the source, the different mass eigenstates accumulate phases 
that depend on their masses, with the differences in these phases governing the beating among the 
mass eigenstates. If the initial state of definite flavor was created as a three-momentum eigenstate (a 
simplification, as one should employ wave packets), the amplitude for finding the state in its original 
flavor at time t and distance ct from the source is 

A{uf ^ Uf) = e^^^ Y: Tf,e-^^^%. (4) 



The time evolution of this state is governed by the Ei = ^JfP + mf, given that a momentum eigenstate 
was assumed, which we expand for small neutrino masses as p + mf/2p. The flavor of the state thus 
evolves if some of the mf differ. The electron neutrino survival probability downstream can be evaluated 
from Eq. 

P{ve i^e) = 1 - sin^ 26113 [cos^ 6*12 sin2(A3iL) + sin^ 6*12 sin^ {/\^2L) \ - cos^ 6*13 sin^ 26112 sm^{/^2iL) (5) 
where L is the distance traveled {fi = c = 1) and 
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Similar probabilities can be evaluated for new flavors appearing downstream, with the probabilities 
summing to 1. The mixing angles Oij and mass differences Sml^ (include their signs) can thus be 
determined from the variation of the neutrino oscillation pattern with distance and beam energy. If 
the case of oscillations among three light neutrinos, only two of the three SrriijS are independent, as 
5mi2 + 6m23 + 6m^i = 0. 

Just as photons traveling through matter propagate according to their index of refraction, neutrinos 
traveling through matter acquire effective masses due to coherent forward scattering off surrounding 
particles [1]. The flavor dependence of this phenomenon - ordinary matter contains electrons, but 
not muons or tauons - can lead to striking effects, enhancing neutrino oscillation probabilities relative 
to those seen in vacuum [2j. The phenomenon is a familiar one, an adiabatic level crossing, and is 
particularly dramatic for small mixing angles. For example, an electron neutrino that might be nearly 
coincident with a light neutrino mass eigenstate in vacuum {Ouip = 0) ~ 0) would be coincident with 
the local heavy mass eigenstate at sufficiently high density (^12 (p — > 00) — )■ tt/2), as the effective mass 
contribution of the matter overcomes the 617121 difference that exists in vacuum. For example, in the 
Sun the electron neutrinos produced in the high-density solar core could be nearly coincident with the 
heavy-mass eigenstates. If the neutrino then propagates adiabatically from high density to the vacuum 
that exists at the solar surface - this requires that the local oscillation length always be small compared 
to the scale height of solar density changes - the neutrino would then emerge into vacuum on the heavy- 
mass trajectory, which would be nearly coincident with, say, the muon neutrino. This phenomenon is 
known as the Mikheyev-Smirnov-Wolfenstein (MSW) effect. The oscillations that occur among three 
neutrino species in the presence of static matter is governed by the equation 




^1 \ / K + K 

rp{23)rj.(13)rj.(12) | Q Q rp(12)trj.(13)trj.(23)t + K 

eJ V v; 



where the Wolfenstein potentials in charge-neutral and unpolarized matter are given by 

VJx) = V2GfNJx) 



for charged-current weak interactions and by 
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(7) 
(8) 

(9) 



for neutral current interactions. In these equations A^^e and Nn are the electron and neutron number 
densities of the medium, respectively. As an overall phase does not contribute to the interferences, we 
can take out a term proportional to the identity by defining Ve = Vc + Vn, V^ = Vr = Vn, and subtracting 
Vfj,l from the second term in Eq. ([T]). Hence, the term Vn drops out of the MSW equations for active 
neutrinos - leaving Vc to generate the electron neutrino effective mass in the solar example given above. 
However, if sterile neutrinos are present, Eq. ([T]) must be expanded to include the additional neutrino 
species, and Vn would then contribute to the potential coupling active and sterile neutrinos. Also note 
that there is a SM loop correction [3] 
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where Np is the proton number density of the medium, is the tau lepton mass, 6w is the Weinberg 
angle, and my/ is the mass of the W-boson. This correction is negligible for the phenomena we consider 
here, but may be important in core-collapse supernovae. Similarly we will ignore interactions between 
neutrinos, which again play an important role in supernovae due to the enormous neutrino densities 
produced by neutrino trapping 
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Often matter-enhanced neutrino oscillations, Eq. ([7|, are treated by retaining just two of the three 
flavors. Indeed if the mixing angle ^13 were zero, Eq. ^7]) can be exactly reduced to a two-dimensional 
problem [5] . As two-flavor evolution is numerically simpler it is still widely employed. The contribution 
of the non-zero value of ^^13 to the electron neutrino survival probability can be obtained using the 
formula [6l [7] 

^3x3(^^6 ^ J^e) = COS^^i3(l -4sin2^13a) P2x2(«^e ^ I^e with iVg COS^ ^13) + siu^ ^13 ( 1 + 4 COS^ ^13a) (11) 

where 

« = ■ ^ ■ V2GFN,ir = 0). (12) 
A32 + A31 

(We are envisioning an application to a case like the Sun, where spherical symmetry can be assumed, 
and thus where densities are a function of the radial coordinate only. Here r = is the solar center.) 
This formula is obtained by expanding the full survival probability as a power series in sin ^13. One can 



calculate the survival probability for the two-flavor case in Eq. (11) by solving the equation 



•i^ /^V^-Asi cos 2^12 ^A2isin2ei2 \ ( (.^^ 

'dtWJ [ |A2isin20i2 + A21 cos 2012/ UxJ ' ^ ^ 

where V = y^GpNe^r) cos^ 613 and ipx = cos 6^23 V^/^ — sin6'23V^r- The MSW level crossing, where the 
adiabatic condition is most severe, is the point where the diagonal terms are equal to zero. 
In some situations the matter basis simplifies the solutions; by making the change of basis 




cos 6'm(r) -sin6'mfr 



sin6'™(r) cos Omir 




(14) 



the flavor-basis Hamiltonian of Eq. (13) can be instantaneously diagonalized. The cosine of this matter 
mixing angle is 

cos2<)„.(.r) = - , (^BE-co.2B,,) ^^^^ 



^{2BE - cos 2^12)2 + sin^ 2^12 
where 

B = V2GfN^ cos^ (16) 

(This result corresponds to the more familiar two-flavor relationship between flavor and local-mass 
eigenstates when 6'i3 — )■ 0.) 

Note that matter effects depend on the sign of Sm'^, i.e., on the neutrino mass hierarchy. Furthermore, 
as the matter background is not CP-symmetric (no antiparticles), matter effects induce an effective CP 
violation. However, in the absence of sterile neutrino mixing, the CP- violating phase 6c p does not 
impact neutrino oscillations even in the presence of a matter background [HI E] . 

3 Measuring 612: Solar and Reactor Neutrinos 

The Sun and about 80% of the visible stars produce their energy by the conversion of hydrogen to 
helium via 

2e- +4p^^He + 2z/e + 26.73 MeV. (17) 

Consequently stars, including our Sun, are prodigious sources of neutrinos. This conversion proceeds 
primarily through the pp chain in lower mass, cooler stars like our Sun, and primarily through the CNO 
cycles in heavier mass, hotter stars. The pp chain, responsible for ~ 99% of solar energy production, is 
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Figure 1: The left frame shows the three principal cycles comprising the pp chain (ppl, ppII, and ppIII), 
the neutrinos sources associated with these cycles, and the theoretical branching percentages defining 
the relative rates of competing reactions, taking from the GS98-SFII SSM Also shown is the minor 
branch ^He+p — >■ ^He+e"'"+z/e, which generates the most energetic solar neutrinos. The right frame 
shows the CN I cycle, which produces the ^^N and ^^O neutrinos. Figure from [TT] . 



comprised of three principal cycles, ppl, ppll, and ppIII, each of which is associated with a distinctive 
neutrino source. The competition among these three cycles in very sensitive to temperature. The CN I 
cycle contribution to solar energy generation depends on pre-existing C and N to catalyze reactions, and 
thus is sensitive to metallicity as well as temperature. Despite its minor role in solar energy generation, 
the CN I cycle produces neutrino fiuxes from the (5 decays of ^^C and ^^N that are potentially measurable. 
The reaction chains are depicted in Fig. [T| including pp-chain branching ratios from a recently updated 
Standard Solar Model (SSM), GS98-SFII [10]. See Ref. [H] for a recent comprehensive review of solar 
neutrino physics, to supplement the summary provided here. 

Pontecorvo first suggested that neutrinos might be detectable by radiochemical means in a large 
volume of chlorine-bearing liquid (see the reproduction of his original paper in Ref. [I2])- This idea 
was subsequently developed by L. Alvarez and later by R. Davis, Jr. In 1964 Davis and his theoretical 
colleague J. Bahcall argued in back-to-back Physical Review Letters [131 E] that a successful solar 
neutrino experiment could be done, given a sufficient detector volume and a deep underground location 
to reduce cosmic-ray-associated backgrounds. They gave as the scientific motivation for their proposal 
"..to see into the interior of a star and thus verify directly the hypothesis of nuclear energy generation..." 
First results from the CI detector were reported in 1968 [15], and it soon became apparent that the 
measured solar neutrino fiux was about one third that predicted by the SSM of the time. 

The source of this discrepancy was unclear for many years because of the large contribution of the 
high-energy neutrinos to the CI capture rate: the fiux of these neutrinos is so sensitive to the solar core 
temperatures that a 5% reduction in that temperature would have brought theory into agreement with 
measurement. The situation began to clarify around 1990, when three new experiments produced data. 
The converted proton-decay water Cerenkov detector Kamiokande II reported the results of 450 days 
of solar neutrino running in 1989, establishing that the ®B fiux was ~ 46% the SSM model prediction 
|16j . in approximate agreement with the Davis result. The signal in this detector was the Cherenkov 
light produced by recoiling electrons after neutrino scattering. This experiment was the first to detect 
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neutrinos in real time and the first to establish their solar origin by correlating the direction of electron 
scattering events with the position of the Sun. First results from SAGE and GALLEX, radiochemical 
detectors employing ''^^Ga as a target because of its low threshold and consequent sensitivity to pp 
solar neutrinos, were reported in 1991 and 1992, respectively. The 1991 SAGE result, based on the 
first five extractions of the daughter nucleus ^^Ge, was an upper bound on the counting rate of 79 
SNU at 90% c.l.[I7] (SNU = solar neutrino unit, 10"^^ captures/target atom/sec). The 1992 GALLEX 
result, based on an initial 14 extractions, was 83 ± 19 (stat) ± 8 (syst) SNU [IB]. Both rates were 
uncomfortably close to the minimum astronomical value of 79 SNU [19], which assumes only a steady- 
state Sun and standard-model weak interaction physics. GALLEX and its successor GNO operated 
until 2003 [201 121]; SAGE remains an active experiment [22]. Because the chemistry of ''^Ge extraction 
from Ga is considerably more complicated than that of ^^Ar extraction from CI, efforts were made to 
calibrate the overall efficiency of both Ga experiments with intense ^^Cr and ^^Ar neutrino sources (~ 
0.5 MCi) [25] [21]. Four such calibration experiments were performed, yielding a weighted average of 
measured to expected "^^Ge of 0.87 ± 0.05 [22]. A counting rate of 66.1 ±3.1 SNU [22| was obtained 
by combining all results from SAGE, GALLEX, and GNO through 2009. This rate can be compared 
to the GS98-SFII SSM prediction of 126.6 ± 4.2 SNU [lO]. 

In the early 1990s, particularly with the uncertainties that accompanied the initial results from 
Kamioka II and the Ga experiments, no individual experiment required a non-astrophysical solution 
of the solar neutrino problem. But in aggregate, the CI, Kamioka II, and Ga experiments indicated a 
pattern of neutrino fluxes that was not compatible with any adjustment of the SSM. This was pointed 
out in a series of papers by various authors as the data became more constraining, eventually establishing 
with a high degree of confidence that the observed rates were inconsistent with the assumptions of a 
steady-state Sun producing undistorted fluxes of solar neutrinos [23 I2SI [23 I2SI 121] ■ 

These early results combined with the spectacular success of the SSM in another context - reproduc- 
ing the sound speed profile deduced from helioseismology to sub-1% accuracy throughout much of the 
solar interior [19] - made a particle-physics solution to the solar neutrino problem more likely. This is 
turn helped motivate a new generation of spectacularly capable active detectors. These experiments are 
Super-Kamiokande, a 50-kiloton water Cherenkov detector from which the ^B flux was determined to 
better than 3%; the Sudbury Neutrino Observatory (SNO), which employed a central 1-kiloton volume 
of heavy water in order to measure three neutrino reaction channels with different sensitivities to flavor; 
and Borexino, which uses liquid scintillator, extending real-time, event-by-event counting to low energy 
fluxes such as the ^Be and pep neutrinos. Both Super-Kamiokande and Borexino continue to operate. 

Results from the SNO detector were definitive in proving new particle physics was responsible for 
the solar neutrino deficit Ray Davis first discovered. Like Super-Kamiokande, SNO was able to detect 
elastic scattering (ES) events 

Z/a; e" Z/a; + C" (18) 

which, due to the charged-current contribution to v^. scattering, has a relative sensitivity to i/gS and 
heavy-flavor neutrinos of ~ 7:1. In addition, SNO could detect z/gS through the charged current (CC) 
reaction 

z^e + d-^p + p + e~, (19) 

by observing the Cherenkov light from the recoiling electron. Unlike the ES case, in the CC reaction the 
e~ carries away most of the energy release, so that the electron recoil spectrum becomes an important 
test of the neutrino spectrum. SNO's third reaction was the neutral current (NC) breakup reaction 

v^{^^)^d^v^{^^)^j)^n. (20) 

This reaction is independent of flavor and thus sensitive to the total flux of neutrinos: by detecting the 
neutron, one can determine a rate that corresponds to an integral over the neutrino spectrum above the 
2.22 MeV deuterium breakup threshold. To exploit this reaction, great care had to be taken to reduce 
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Figure 2: The SNO CC, ES, and NC and Super-Kamiokande ES constraints on the z/g and u^r compo- 
nents of the flux are shown. The SSM total flux is also shown. All results are all consistent with a 
terrestrial flux of solar neutrinos that is about one-third i/gS and two-thirds z/^r |31]. Reproduced with 
permission from [32] . 



background sources of neutrons from natural radioactivity and from cosmic ray muons. SNO's great 
depth, two kilometers below the surface, was helpful in addressing the latter. The SNO experiment 
was carried out in three phases in which distinct methods were used in the NC channel, neutron 
capture on deuterium, on CI introduced as salt in the heavy water, and in an array of '^He proportional 
counters constructed in the central detector. The SNO results for the CC, ES, and NC channels, the 
Super-Kamiokande ES results, and the predictions of the SSM circa 2005 are shown in Fig. |2} The 
data are consistent, showing that approximately two-thirds of the solar neutrinos arriving at earth are 
heavy-flavor neutrinos |30l[3T]. The total flux is consistent with the SSM, given assigned uncertainties. 

Table [T] summarizes the properties of the solar neutrino experiments carried out to date. 

The solar neutrino demonstration of neutrino oscillations was confirmed in the KamLAND reactor 
neutrino experiment [33] . The KamLAND experiment used neutrinos from 55 Japanese nuclear reactors 
with a flux-weighted average baseline of 180 km, a distance fortuitously appropriate for observing 
oscillations governed by the solar 5m\^. The KamLAND results provided a precise determination of 
(5m^2 ^ ('''•^l-aig) X 10~^ eV^, eliminating an ambiguity that existed in analyses that included only 
solar neutrino data: both the LMA (large mixing angle) solution corresponding to 5mi2 ~ 8 x 10~^ eV^ 
and the LOW solution, 5m\2 ~ 10^'^ eV^, were allowed at 67% c.l., prior to KamLAND measurements. 
KamLAND results eliminated the LOW solution. 

Neutrino propagation in the Sun is influenced by matter, which alters the relationship between the 
local mass eigenstates, effectively rotating the local mixing angle toward tt/2 with increasing density. 
The electron neutrino survival probability at Earth, averaged over the solar core production region 
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Table 1: Solar neutrino experiments. 



Experiment 

Sr 


Type 


Method 


Dates 


Threshold(MeV) 


Main Fluxes 




i. CAj'^tl ^1 J-Vj J. J. J. 


CC 


1968-2002 


0.811 


^B, ^Be 


TC PI m i olcpi n H p 

J. i.CtjlJ.AlWl.VCtiJ.J.'—tVj 


r\ ft i VP 


ES 


1987-1995 


7.0-9.0 


^B 


SAGE 


rriHiofVipm 

i. Ctj*— 4.1 ^1 J-Vj J. J. J. 


CC 


1 QPO-nrp^pnt, 


0.233 


nn '''Bp 


GALLEX/GNO 


radiorhem 


CC 


1991-2003 


0.233 


nn ^Bp 


Super-Kamiokande 


active 


ES 


1996-present 


4.0-7.0 


8B 


SNO 


active 


ES/CC/NC 


1999-2006 


4.0-7.25; 2.22(NC) 


«B 


Borexino 


active 


ES 


2008-present 


~ 0.8 


^Be, pep, ^B 



(which removes interference terms between the local heavy and light mass eigenstates), is 

P2x2{^e ~^ = ^ + ^ cos 26*12 (cos 26'i)so„^ce (1 " '^Phop) ■ (21) 

Here (cos 26*4) source is the local mixing angle at the point of neutrino production, suitably averaged over 
the neutrino-producing region of the solar core; the superscript D indicates that the neutrino has been 
observed during the day, so that matter effects associated with passage through the earth play no role; 
and the subscript 2x2 indicates that subdominant oscillations involving the third flavor have been 
neglected. Phop is the probability that the neutrino hops from the heavy-mass trajectory to the light- 
mass trajectory on traversing the level-crossing point [3^ |35]. This accounts for possible non-adiabatic 
behavior: if the density in this region is treated linearly, one finds 



hop 



e-">/2 (22) 



where 7c is proportional to the ratio of the Sun's density scale height to the local oscillation length at 
the crossing point. (Thus large 7c corresponds to small changes in the solar density over an oscillation 
length, the adiabatic limit where Phop — J- 0.) 

Solar neutrinos detected at night go through the Earth, hence they would feel additional matter 
effects [36]. The night-time solar electron neutrino survival probability, averaged over the Earth-Sun 
distance, is given by EH] 

P2>c2{^e ^ Z/e) = ^ + \{cOs2e^) source (1 - 2Phop) COs2de COs2(^e " ^12) (23) 



where 9e is the matter mixing angle inside the Earth. In writing Eq. (|23|), matter density in the mantle 

(24) 



of the Earth is assumed to be constant. The day-night asymmetry A is defined as 



2 pN j^pD- 

The terrestrial matter effects on solar neutrinos that can generate day/night differences are quite small: 
based on current neutrino parameters obtained from global analyses, a ~ 3% day-night effect is expected 
in ES measurements of the ®B spectrum. The latest available Super-Kamiokande results remain consis- 
tent with no effect at 2.6 a [39j, though it is quite possible that with the data now being accumulated 
in Super-Kamiokande IV low-threshold running, this conclusion may change. The corresponding SNO 
result is less restrictive, consistent with the null hypothesis of no day/night effects at 61% c.l. [31] . 
The predicted day/night effect in Borexino is tiny, consistent with the collaboration's determination A 
= 0.001 ± 0.012 (stat) ± 0.007 (syst) [30]. (The absence of day/night effects in Borexino's ^Be flux 
measurement is an alternative means of eliminating the LOW solution.) 
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In contrast, matter effects associated with passage through the Sun are large, with the additional 
suppression being substantial for the higher energy neutrinos. Because matter effects either increase 
of decrease the local splitting between the mass eigenstates, depending on the hierarchy, solar matter 



effects allow one to establish that Sm2i is positive (m2 > rrii). Consequently B of Eq. (16) is positive. 
When one uses neutrino oscillation parameters taken from current global analyses, one finds that low- 
energy solar neutrino fluxes (pp, pep, ''Be) reside in a region where vacuum effects are dominant, 

pvaccum ^1-1 giu^ 2^12 ~ 0.57 (25) 

while oscillation probabilities for the high-energy portion of the ^B neutrino spectrum measured by 
SNO and Super-Kamiokande are substantially affected by matter 

phigh density ^ ^-^2 _ q ^26) 

These expectations are clearly apparent in the data: the pp flux determined in luminosity-constrained 
global analyses and the ^Be and pep fluxes determined by Borexino [HI all show the expected 
turn-up in P,,^ at low E^. (See Fig. |3}) 

One of the goals of both SNO (CC/ES) and Super-Kamiokande (ES) has been to see the matter 
effects imprinted on the scattered electrons. SNO's low-energy-threshhold (LET) analysis [31] and 
Super-Kamiokande IV running both employ total-energy thresholds of ~ 4.0 MeV. This provides sen- 
sitivity to a region of the spectrum where the matter-to- vacuum transition might become apparent. 
The Pjy^s deduced from the SNO LET analysis are included in Fig. |3} While remaining consistent with 
global analysis predictions, the best flt is trending away from expectations. 

The three-flavor analysis based on all solar data supplemented by KamLAND yields 



sin^ 612 = 0.308 ± 0.014 6ml^ ~ (TAltoH) x 10~ W sin^ ^13 = 0.025tomt (27) 

These results are in excellent agreement with the global analysis that takes into account all neutrino 
data, including beam experiments and new reactor experiments. The Bari |13] and Valencia [H] analyses 
yield, respectively, 

. 0.3061°:°}^ r (7.54^°:i) . 0.0243t°:°°r, 

sin^ 9^2 = \ 5ml-^ ~ <^ x 10-^eV^ sin^ = \ (28) 

i 0.3191H1^ ^ (7.62i°i^) i 0.0248i°:°°I« 

The solar neutrino analysis supplemented by KamLAND produces a best-value estimate of 6*13 that is 
in very good agreement with the Bari and Valencia results. However, the global analyses provide a 
more certain determination of 613: this reflects the impact of the recent Daya Bay, RENO, and Double 
Chooz reactor experiments, which tightened the error bars on ^13 by almost an order of magnitude. We 
discuss these experiments in Sec. [5j 

The original goal of solar neutrino spectroscopy was to use neutrinos to test the hypothesis that 
the Sun is powered by the pp chain and to measure the solar core temperature to high accuracy. This 
program was delayed by the long effort to identify the origin of the solar neutrino problem and then 
to accurately determine the parameters governing neutrino oscillations. Consequently today, given 
the few-percent determinations made of the relevant mixing angles and 5m\2i it has flnally become 
possible to use neutrinos as quantitative astrophysical probes of the Sun and the SSM. In addition to 
weak interaction uncertainties, the SSM depends on the accuracy with which an additional ~ 20 input 
parameters are known. These parameters are associated with SSM's nuclear and atomic physics and 
with properties such as the Sun's initial metallicity, mass, age, luminosity, and the diffusion coefficient. 
Decades of laboratory and theoretical work on the nuclear and atomic physics input has helped make 
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the SSM a predictive model, where output uncertainties are known. The laboratory astrophysics of 
the SSM was advanced by community efforts like that leading to the analysis of cross sections and 
S- factors published in 1998 [17] (Solar Fusion I). The key uncertainties identified in that analysis led 
to a decade of further work on critical cross sections, such as that for ^Be(p,7), and to the recent re- 
evaluation of cross sections summarized in [18] (Solar Fusion II). Current SSM predictions of neutrino 
and helioseismic observables have uncertainties based on evaluated input laboratory data. Discrepancies 
in SSM predictions that significantly exceed model uncertainties could indicate problems with the 
underlying astrophysics assumptions of the SSM. 

The key discrepancy in current SSM predictions is associated with its initial metallicity (abundances 
of elements heavier than hydrogen and helium). The overall solar metallicity scale is usually determined 
from analyses of solar photo-absorption lines. Recently more realistic models of the photosphere were 
constructed to account for its dimensionality and granularity, effects ignored in traditional ID slab 
models. When applied to the photospheric data, these new models yield an abundance ratio of metals 
to hydrogen (Z/X)0 = 0.0178 [19], compared to the previously accepted value of (Z/X)© = 0.0229 
|50j . Recent SSM calculations have been carried out for both metallicities [1^, while also incorporating 
the new results from Solar Fusion II. The helioseismic predictions of the low-metallicity SSM are in sig- 
nificant disagreement with observation. The neutrino flux predictions of both models are in reasonable 
agreement with the data, though it is apparent that a solar core metallicity set to the average of the 
low- and high-metallicity models would improve the agreement. 

This has opened up two directions of study. First, one of the central assumptions of the SSM is 
that the Sun was homogeneous at the start of the main sequence. There is little observational support 
for this assumption, given that it is known, late in the evolution of the solar disk, that the process of 
planetary formation scoured 50-90 of excess metal from the disk. New solar models have been built 
to take account of the possible accretion of metal-depleted gas onto the proto-Sun, late in its evolution 
when such accretion would alter convective-zone abundances but not affect the core. It has been found 
that both solar neutrino and helioseismic data severely constrain such accreting models [10] . Second, 
the volatile elements CNO play a key role in the solar abundance problem. The SSM produces about 
1% of its energy through the CN cycle, yielding significant fluxes of ^^N and ^^O solar neutrinos. It has 
been shown, because of the CN cycle's additional linear dependence on the abundance of C+N (these 
metals catalyze the proton burning), that the C+N abundance of the solar core could be determined 
to an accuracy of about 10% from the CN neutrino fluxes [51]. This would be very significant given 
the current 30+% discrepancy between the low- and high-metallicity SSMs. The CN neutrinos can be 
measured in very deep scintillation detectors such as SNO+, and possibly also in shallower detectors 
such as Borexino, if strategies for suppressing cosmic-ray muon backgrounds can be further developed. 

There are other gaps in our program of solar neutrino spectroscopy. While the pp flux is tightly 
constrained by the solar luminosity constraint, given existing ''Be and ^B neutrino flux measurements, 
this constraint itself should be checked: are the Sun's neutrino and photon lumninosities consistent, as 
predicted in steady-state solar models? Recently the Borexino collaboration succeeded in detecting the 
pep neutrinos [12], though the accuracy of the measurement is so far limited. This flux should track 
that of the low-energy pp neutrinos, the dominant solar neutrino flux arising directly from p+p fusion. 
While the pp flux is constrained by the results from the Ga experiments, no real-time detection has 
been made. The relation between the photon and neutrino luminosities in steady-state models is [52] 

^, = T.iQ-{m<p. (29) 

where d is the average Sun-Earth distance (1 a.u.), Q is the energy released in the nuclear reactions, {E) 
is the average energy loss by neutrinos, and (pi is the neutrino flux at earth, coming from the reaction of 
type i. (There are corrections to this equation due to the necessary inequalities between different fluxes 
to take into account the order in which nuclear reactions take place [53t [29]). A test of this relation 
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would require a precise measurement of both the pp- and CNO-chain neutrino fluxes. Such a test would 
place stringent limits on new neutrino species, while also checking the assumption of steady-state solar 
burning. 



4 Measuring 623: Atmospheric and Accelerator Neutrinos 

The hadronic cosmic ray flux incident on Earth is dominated by energetic protons with fewer numbers of 
heavy nuclei. When the primary cosmic rays collide with the oxygen and nitrogen in Earth's atmosphere 
they produce secondary pions, kaons, and muons. Atmospheric neutrinos come from the decay of these 
secondaries. For energies less than 1 GeV all these secondaries decay in the atmosphere 

/i^ ^ + Ue{Ve)+lyM- (30) 

Consequently one expects the ratio 



r = + v^l (z/^ + z/^) (31) 

to be approximately 0.5 in this energy range. This relationship can be checked as a function of distance 
up to ~ 13000 km, by observing atmospheric neutrinos in a deep detector as a function of zenith angle. 
It has been customary to use a ratio-of-ratios 

^ ^ [(^e + t^e)/(t^M+^A')]data (32) 
[(z/e + Fe)/K+l^M)]Monte Carlo 

in this check, normalization the measured r to a Monte Carlo calculation of the ratio assuming no 
oscillations but including corrections such as geomagnetic effects on the cosmic ray flux. 

Early atmospheric neutrino experiments were inconsistent with one another, with some finding i? ~ 1 
while others finding diminished values. Significant indications of an anomaly came from the 1MB [5lj 
and Kamiokande [55] detectors. 1MB first noticed a possible deficit of neutrino-induced muon events in 
1986, while Kamiokande established a deficit in excess of 4cr by 1988. By 1998 the anomaly was also 
apparent in data taken from the Soudan detector and from Super-Kamiokande 

In early experiments it was difficult to adequately explore R because the lack of statistics prevented a 
detailed analysis as a function of zenith angle. This changed with the construction of Super-Kamiokande, 
which provided a fiducial volume of ~ 20 ktons. In an initial analysis the Super-Kamiokande Collabo- 
ration found 

^ _ f 0.61 ± 0.03(stat) ± 0.05(sys) sub — GeV events, fully contained 

y 0.66 ± 0.05(stat) ± 0.08(sys) multi — GeV events, fully or partially contained 

The 1998 Super-Kamiokande analysis, based on 33 kton-years of data, showed a zenith-angle dependence 
for R inconsistent with theoretical calculations of the atmospheric fiux, in the absence of oscillations. 
The distance-dependent muon deficit was consistent with oscillations governed by neutrino mass dif- 
ferences within the range 5 x 10-^eV^ < \bm\^\ < 6 x lO'^eV^ [SB]. This differed from that later 
determined for solar neutrinos. The collaboration concluded that the measurements were consistent 
with the two-fiavor oscillation z/^ — v^^ providing evidence for massive neutrinos [56] . 

Super-Kamiokande I collected ~ 15,000 atmospheric neutrino events over approximately five years 
of running. The collaboration's zenith-angle analysis of the data placed the first oscillation minimum at 
L/E 500 km/GeV, so that Lq ~ 1000 km for a 1 GeV muon neutrino. Results based on 3903 days of 
data from Super-Kamiokande I-IV, including 1097 days from the current phase IV, were incorporated 



12 



Table 2: Values of sin^ 26'i3 measured at various experiments. The value quoted for the 
MINOS experiment corresponds to 2 sin^ ^23 sin^ 2^13. 

Experiment Measured sin^20i3 



Daya Bay 0.089 ± O.OlO(stat) ± 0.005(syst) 

RENO [70] 0.113 ± 0.013(stat) ± 0.019(syst) 

Double Chooz [HH] 0.109 ± 0.030(stat) ± 0.025(syst) 

T2K (assuming ^23 = vr/4) [72] 0.104lg:[]^g 
MINOS [73] (see note) +0.04l;g:g^I (normal) 0.079;g:g^;^ (inverted) 



5m|2 = 


(2.66l°i^) 


X lO-^eV^ 


(la) 


sin^ 6'23 = 




U.^ZO_o.o34 


(90%c.l.) 


6ml^ = 


(2.66l°:M) 


X lO-^eV^ 


(la) 


sin^ 6*23 = 




n p;7p;+o.055 


(90%c.l.) 



in a three-flavor analysis in which sin^ ^13 was set to 0.025, based on Daya Bay, RENO, and Double 
Chooz results. The analysis yielded [STJ |58] 

Normal hierarchy : 
Inverted hierarchy : 

The significant difference between the normal- and inverted-hierarchy best values for sin^ ^23 reflects 
the flatness of the fit around sin^ ^^23 ~ 0.5, with shallow local minima appearing above and below 
this value. 

Accelerator experiments were conducted by the K2K [59] and MINOS [6OI [61] collaborations, with 
baselines designed to be sensitive to the scale probed with atmospheric neutrinos. The MI- 

NOS beam-neutrino analysis, which includes atmospheric neutrino data from that detector, determined 
ISml^l = 2.39t 

aio [62], a value consistent with the somewhat broader range allowed by the at- 
mospheric results. The sign of not been determined, though one can see from the results 
above that there is some sensitivity to the hierarchy in atmospheric neutrino fits. There are several 
pending proposals to determine the hierarchy: new, long-baseline neutrino beam experiments (some 
of these depend on joint analyses of beam and atmospheric neutrino data), high-density arrays to en- 
hance the sensitivity of large detectors like IceCube to atmospheric neutrinos, and new reactor-neutrino 
experiments with very large detectors to increase sensitivity to subdominant oscillations. 



5 Measuring 6i^: Reactor and Accelerator Neutrinos 

Initial comparisons of the solar neutrino and KamLAND reactor neutrino experiments revealed hints of 
a non-zero value of ^13 [631 El]- The solar neutrino+KamLAND analysis performed by the KamLAND 
Collaboration in 2010 yielded sin^ ^13 = 0.020t[]:J^ [65] . with a nonzero value suggested at 79% c.l. 
Indications of electron neutrino appearances from accelerator-produced off-axis muon neutrinos in the 
T2K [2S] and MINOS [ST] experiments also suggested relatively large values of ^13. In 2012 the Double 
Chooz reactor collaboration announced data that suggested a nonzero value for sin^ 26*13 at 94.6% c.l. 
|68j . followed in quick succession by definitive results from the Daya Bay [69] and RENO collaborations 
[70] . The various determinations of sin^ 2^13 are summarized in Table [2} 

Reactor experiments are disappearance experiments, with ^13 determined from the magnitude of the 
electron antineutrino flux loss with distance. The disappearance probability, taken from Eq. (|5]), is 
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cos^ 6*12 sin^(A3iL) + sin^ 6*12 sin^(A32L) +cos'' 6*13 sin^ 26*12 sin^(A2iL). (34) 



If the experiment is done close to the reactor, A21L ~ so that the "solar" oscillation represented by 
the second term can be ignored. Using the fact A32 ~ A31, suggested by the solar neutrino experiments. 
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one obtains the disappearance probability 



1 - P(u^ Z/p 



sin^2^i3 sin2(A3iL). 



(35) 



Consequently the very short-baseline reactor neutrino experiments, such as Daya Bay, Double Chooz, 
and RENO, unambiguously measure ^13, requiring no knowledge of other neutrino mixing angles or the 
mass hierarchy. Furthermore, the multiple detector configurations currently employed in the Daya Bay 
and RENO experiments also minimize the effects of reactor neutrino spectrum uncertainties. 

The situation is rather different for longer-baseline experiments searching for the appearance of 
electron neutrinos in a flux of muon neutrinos. For longer baselines neutrinos travel trough the Earth, 
hence matter effects need to be included. The resulting rather complicated expressions, which depend 
on many of the neutrino parameters, can be calculated in a series expansion [71]. The appropriate 
appearance probability in the lowest order is given by 



sin^ 29i3 sin^ 623 



sm 



A 



31 



V2 



+ 0{g), 



(36) 



where 
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5im} 
6m 



21 



0.03. 
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(37) 



The next order correction in g to the electron neutrino appearance probability brings in an additional 
dependence on the CP-violating phase in the neutrino mixing matrix. The denominator of the term 



multiplying the oscillating term in Eq. (36) depends on the sign of Sm , i.e. the mass hierarchy of the 



neutrinos. Consequently, appearance experiments such as T2K [72] and MINOS [73] cannot disentangle 
^13 from other observables such as the mass hierarchy, 623, and the CP- violating phase. Note that an 
intermediate-baseline (L ~ 60 km) reactor antineutrino experiment will also have some sensitivity to 
the neutrino mass hierarchy [75j . 



6 Final Remarks 

The puzzle that arose from early solar neutrino measurements had an interesting resolution, still incom- 
plete. We have reviewed what has been learned in the past 15 years about the properties of neutrinos 
from oscillations in vacuum and in matter, using astrophysical neutrino sources and terrestrial beams. 
A great deal of progress has been made, with three mixing angles and two mass differences now ac- 
curately determined, and with the sign of 5m^2 deduced the effects of solar matter on solar neutrino 
oscillations. Yet there are several important properties of the light neutrinos that we have not yet been 
able to determine: the overall scale of the neutrino masses, the mass hierarchy, the Dirac CP-violating 
phase in the mixing matrix, and neutrino charge-conjugation properties. We have the tools to make 
further progress. Long-baseline neutrino beam oscillation experiments could distinguish between the 
normal and inverted mass patterns and determine the CP phase. As we have noted, there may be alter- 
native strategies in which reactor and atmospheric neutrino sources are exploited to probe the neutrino 
mass hierarchy. A variety of cosmological tests are sensitive to the absolute scale of neutrino mass: the 
distinctive imprint of neutrino mass on the growth of large-scale structure, affecting both the scale and 
red-shift dependence of that growth, produces changes in the power spectrum that are roughly an order 
of magnitude larger than one might naively guess, based on the contribution of neutrinos to the mass 
density. There are a variety of ongoing projects to improve laboratory sensitivity to neutrino mass, 
such as the high-sensitivity spectrometer KATRIN for measuring tritium beta decay and the MARE 
micro-calorimeters for measuring ^^^Re beta decay, as well as new technologies for overcoming resolu- 
tion challenges, such as the Project 8 extension to KATRIN. Finally, the exciting question of whether 
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Nature exploits the neutrino's unique Majorana mass mechanism is stimulating next-generation efforts 
in double beta decay. As progress is made on all of these fronts, neutrino oscillations will remain one 
of our most important and most versatile tools from probing the new properties of neutrinos, and thus 
for seeking new physics beyond the SM. 
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under ,de-sc00046548 (UC Berkeley) and DE-AC02-98CH10886 (LBL), and by the UC Office of the 
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